Relativistic theory of the double photoionization of helium-like atoms 
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A fully relativistic calculation of the double photoionization of helium-like atoms is presented. The 
approach is based on the partial-wave representation of the Dirac continuum states and accounts 
for the retardation in the electron-electron interaction as well as the higher-order multipoles of 
the absorbed photon. The electron-electron interaction is taken into account to the leading order 
of perturbation theory. The relativistic effects are shown to become prominent already for the 
medium- Z ions, changing the shape and the asymptotic behaviour of the photon energy dependence 
of the ratio of the double-to-single photoionization cross section. 
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I. INTRODUCTION 



The ejection of two electrons caused by absorbtion 
of a single photon is one of the fundamental few-body 
processes in atomic physics. The process is called dou- 
ble photoionization or, less frequently, photo-double- 
ionization. The characteristic feature of this process is 
that it proceeds exclusively through the electron-electron 
interaction. Because of this, double photoionization has 
long been used as a testing ground for understanding of 
the electron correlation phenomena. 

The traditional system for studying double photoion- 
ization is the helium atom, for which numerous exper- 
imental and theoretical investigations have been per- 
formed during the last four decades (for a recent re- 
view, see Ref. [l|). Most widely studied is the ra- 
tio of the double-to-single photoionization cross sections 
R = a ++ 1 a + as a function of the energy of the incoming 
photon uj. Early calculations of this ratio were mainly 
concentrated either at the near-threshold region of the 
photon energies {uj > lu ci 79 eV for helium, with w cr 
being the double ionization energy), where Wannier the- 
ory is shown to be applicable 0, H| , or at the asymptotical 
nonrelativistic region uj ct <§C w <C to 0-11 (where m is the 
electron mass). The intermediate region of the photon 
energies (uj ~ 200 eV for helium, where R(uj) reaches its 
maximum) turned out to be much more difficult for an 
accurate theoretical treatment. Reliable theoretical pre- 
dictions in this region were obtained only in the end of 
1990th by means of sophisticated many-body techniques 
(notably, the close-coupling methods) |7HTu|. 

Since the helium atom is an essentially nonrelativistic 
system, any relativistic effects in its double photoion- 
ization are considered to be of little importance at the 
present level of experimental precision. However, the re- 
cent experiments on the double Af-shell pho toionization 
of moderately heavy atoms up to silver jll| - [l5| demon- 
strated significant enhancements (of about factor of five 
for silver) of the cross section as compared with results 
of nonrelativistic calculations. A natural candidate for 



explaining this enhancement would be relativity, which 
obviously cannot be disregarded when dealing with the 
deeply bound electron states in silver. Because of the 
relative isolation of the A'-shell electrons from the outer 
electrons, the double if-shell ionization of a heavy atom 
is often compared to the double photoionization of the 
corresponding helium-like ion. However, an accurate the- 
oretical treatment of double A"-shell ionization should in- 
clude both the relativistic effects on the inner-shell elec- 
trons and the electron-correlation effects induced by the 
outer shells. Such a calculation is rather difficult and has 
not been performed so far. 

In contrast to many-electron atoms, the helium-like 
ion is a relatively simple system for which an ab ini- 
tio description is feasible. So far, there have not been 
any direct measurements of the double photoionization 
of helium-like ions. However, with the advent of new 
powerful light sources, such as the free-electron laser 
(FLASH) in Hamburg, the Linear Coherent Light Source 
(LCLS) at Stanford, and the X-ray Free Electron Laser 
(XFEL) at Hamburg, the experimental study of various 
photoabsorbtion processes is going to be possible for a 
great variety of ions in different charge states [l6| . Mea- 
surements of the double photoionization in moderately 
heavy helium-like ions would allow us to effectively test 
our understanding of the relativistic electron-correlation 
effects. 

An alternative approach to the investigation of the 
double photoionization of highly charged ions is to study 
it in the inverse kinematics, through the radiative dou- 
ble electron capture. For bare oxygen, such measurement 
was recently accomplished in Ref. [17j . The experimen- 
tal upper bound on this process for bare uranium was 
reported in Ref. fTsj j . 

The goal of the present investigation is to perform an 
ab initio relativistic calculation of the double photoion- 
ization of a helium-like atom. The electron-electron cor- 
relation (both on the initial and the final states) will be 
taken into account to the leading order of the perturba- 
tion theory. To the given order of perturbation theory, 
the treatment is exact, i.e., includes all multipoles of the 
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absorbed photon, the retardation in the electron-electron 
interaction as well as the interaction of the electrons 
with the nucleus without any expansion in the binding 
field. The treatment is gauge invariant, both with re- 
spect to the gauge used for the absorbed photon and the 
gauge of the electron-electron interaction. The higher- 
order electron-correlation effects omitted are estimated 
by comparing the present numerical results with the ex- 
perimental data available for helium. 

The remaining paper is organized as follows. In Sec. [TT] 
we present a short summary of the rclativistic formulas 
for the single photoionization, which forms the basis for 
our treatment of the double photoionization. In Sec. IIIII 
we describe the QED theory of the double photoioniza- 
tion to the lowest relevant order of the perturbation the- 
ory. Sec. IIVI presents details of the calculation. Numeri- 
cal results are presented and discussed in Sec. El 

The relativistic units (K = c = m = 1) are used 
throughout this paper. 



II. SINGLE PHOTOIONIZATION 

In this section we summarize the relativistic formulas 
for the single photoionization, as they build the basis for 
the description of the double photoionization. 

Differential cross section of the photoionization of a 
one-electron atom is 

d(7+ = ^LJl | T +| 2 £( e . +u-e)dp, (1) 

where a is the fine-structure constant, u) is the energy 
of the absorbed photon, £j is the initial (bound-state) 
energy, e and p are the energy and the momentum of 
the emitted electron, respectively, e = \J m 2 
amplitude of the process T~t is given by 
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where \niiJLi) denotes the initial Dirac bound state with 
the relativistic angular quantum number Kj and its pro- 
jection /ij, and \prnj is the wave function of the emitted 
electron. The general relativistic expression for the pho- 
ton absorption operator R\ is 



R\ = a. ■ u\ e l 



G(a-k-l) 



(3) 



where a is a three-component vector of the Dirac matri- 
ces, u\ is the polarization vector of the absorbed pho- 
ton, k is the photon momentum, k = k/\k\, and G is the 
gauge parameter. In our treatment, all electron states 
are the eigenfunctions of the Dirac Hamiltonian, so the 
gauge-dependent part of R\ vanishes identically. How- 
ever, we keep the general gauge of the absorbed photon 
in actual calculations in order to check the numerical 
procedure. 

The wave function |pm) = |pm)_ in Eq. (|2|) is the 
one-electron continuum Dirac state with the asymptotic 



momentum p, helicity m = ±!/2, and the "— " asymptotic 
behaviour (i.e., the plane wave modified by the Coulomb 
logarithmic phase plus the incoming spherical wave) [l9| 



\ pm )_ = e <b» 
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as \r — z\ — > co, where rj = Zas/p is the Sommerfeld 
parameter. The free-electron 4-spinors u(p,m) (2(J are 
normalized by the condition uu = u + ^/qu = 1. The ex- 
plicit expression of the electron wave function is 
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where A; are the scattering phases [21|, j 



I = \k+ V 2 | — V 2 7 and \enn) are Dirac continuum states 
with given relativistic angular momentum n and the mo- 
mentum projection /i, normalized on the energy scale 
19]. 

Taking into account that dp — pededQ, the sin- 
gle differential cross section of the photoionization of a 
hydrogen-like atom is written as 
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where the energy of the emitted electron is fixed by the 
energy conservation, e = Ei + lu. 

In our analysis of the double photoionization, we will 
need the cross section of the single photoionization of 
a helium-like atom. In the independent-particle approxi- 
mation, the wave function of the initial two-electron state 



is 



\J M ) = NY, C£Z\ 



X — (\K a fl a )\K b fl b ) - \KbfJ-b)\KafJ-a)) , (J) 

where N = 1/V2 for the equivalent electrons and N = 1 
otherwise. Averaging the cross section over the momen- 
tum projection of the initial state Mo, employing the ex- 
plicit expression for the continuum Dirac state ([5]), and 
integrating over the angles over the emitted electron, we 
obtain the total photoionization cross section of a helium- 
like atom as 
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where e\ = e a + w and Ei = e b + w and we assumed that 
the energy of the photon is sufficient to ionize any of the 
two initial-state electrons. If this is not the case, only 
one of the two terms in the brackets should be retained. 



3 



Pi 




a 



P2 b p 2 b 

FIG. 1: Feynman diagrams representing the double photoion- 
ization of a helium-like atom, a and b denote the bound elec- 
tron states, pi and P2 are the continuum electron states, k 
denotes the incoming photon. Double lines denote electrons 
propagating in the binding nuclear field. 

The electron correlation effects on er + omitted in the 
independent-particle approximation are well studied in 
the literature. In particular, the leading nonrelativistic 
term of the 1/Z expansion of the high-energy asymptotics 
of cr + was obtained in Ref. (22j. However, since we are 
presently interested in the ratio of the double-to-single 
photoionization cross sections a ++ /a + , we prefer to treat 
both these cross sections on the same footing, i.e., to the 
leading nonvanishing order of perturbation theory with 
respect to the electron-electron interaction. 



the double photoionization by many authors, notably in 
Refs. d, Although the formal parameter of the per- 
turbative expansion is 1/Z, the leading-term approxima- 
tion was shown to adequately describe the double pho- 
toionization even for helium. Our present treatment is 
aimed primarily at the helium-like ions, for which the 
perturbative expansion converges much faster than for 
helium. 

The triple differential cross section of the double pho- 
toionization obtained from Eq. ([9]) is 

d 3 a++ An 2 a 
dilidil 2 d£i uj 

x ^2 \ T \ + (Pi m i,P2m 2 ; uj, JoM )\ 2 , 

mi??i2 

(ii) 

where the energy of the second electron is fixed by the 
energy conservation, £2 = £a + £b + uj — E\. The energy 
distribution of the emitted electrons can be conveniently 
parameterized by the fractional energy sharing parameter 

x, 

£\ = m + x (uj — uj cr ) , (12) 
e 2 = m + (1 - x) (uj — u) CI ) , (13) 



III. DOUBLE PHOTOIONIZATION: GENERAL 
FORMULAS 

According to the general rules of quantum field theory 
[23j , the differential cross section of the double photoion- 
ization is 

At: 2 a 2 

da ++ = T++ S(£i + uj - £1 - £ 2 ) dpi dp 2 , (9) 

ui ' 

where (ei,pi) and (£2,^2) are the energy and the mo- 
mentum of the two final-state electrons, respectively, £j 
is the energy of the initial bound state, and is the 
amplitude of the process. To the leading order of per- 
turbation theory, the amplitude of the process is repre- 
sented by the two Feynman diagrams shown in Fig. [T] 
where the antisymmetrization of the initial and the final 
two-electron states is assumed. 

The initial-state wave function is the same as for the 
single photoionization [see Eq. ([7])] and the initial-state 
energy is £j = e a + £b- The wave function of the final 
state is 

\pim 1 ,p 2 in 2 )_ = (\pi7ni) _\p 2 m 2 ) _ 



-|P2m 2 )_|pimi) 



(10) 



Note that in the present approach, we use the wave 
functions that do not include any electron-correlation ef- 
fects (apart from the antisymmetrization); the electron- 
electron interaction enters explicitly into the amplitude 
of the process [see Eq. (|16l) below] . In the nonrelativis- 
tic case, this approach has been successfully applied to 



where uj cv is the threshold value of the photon energy 
(the double ionization energy). In our approximation, 
ui CI = 2m - e a - Eh- 

Let us now consider the single differential cross section 
da ++ 1 ' de\. Substituting the explicit expression for the 
continuum wave function ([5]) into Eq. (jlll) . integrating 
over the angles, and summing over mi and m 2 , we arrive 
at 



At: 2 a 



de\ 



E 
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(14) 



where the amplitude contains only the spherical-wave 
Dirac continuum states. 

The total cross section is obtained as an integral of the 
single differential cross section over a half of the energy 
sharing interval 



m-\-(uj — ~\ 



de\ 



dsi 
da++ 



(uj - u) CT ) / dx , 
Jo dei 



(15) 



The other half of the energy interval corresponds to in- 
terchanging the first and the second electron, which is 
already accounted for by the wave function. 

The general expression for the transition amplitude of 
the double photoionization process can be obtained by 
the two-time Green's function method 231 ] . 
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(16) 



The first term in the brackets corresponds to the electron- 
electron interaction modifying the final-state electron 
wave function (the left graph in Fig. [T|) and the second 
one, the initial-state electron wave function (the right 
graph in Fig. [T]). The summation over P and Q corre- 
sponds to the permutation of the initial and final elec- 
trons, Pe\Pei — {e\z<x) or (£ 2 £i), QaQb — (ab) or (ba), 
and (— l) p and (—1)^ are the permutations sign. The 
summation over n in Eq. (|16p runs over the complete 
Dirac spectrum, A a b = e a — £&, and 1(A) is the electron- 
electron interaction operator, 



= e 2 a^a v D^(oj,x 12 ), 



(17) 



where D 111 ' is the photon propagator, = (l,ct) are 
the Dirac matrices. In the Feynman gauge, the electron- 
electron interaction takes the form 



/ Fcyl, (w)=a(l-a r a 2 ) 



Xl2 



(18) 



where X\% — \xi — xa|. In the Coulomb gauge, the 
electron-electron interaction acquires an additional term, 
which can be expressed as 



I Cou \uj) =I Fcyn (uj) 

(c*i • Vi)(a 2 • V 2 ) 



or 



X12 



(19) 



IV. NUMERICAL CALCULATION 



The final formulas for the single and double photoion- 
ization cross sections are presented in Appendices O and 
151 The calculation of the single photoionization cross sec- 
tion is straightforward, in contrast to that of the double 
photoionization. The major difficulty in the numerical 
evaluation of the double photoionization is the summa- 
tion over the complete spectrum of the Dirac equation. 
In this work we use the approach based on the analyti- 
cal representation of the Dirac Coulomb Green function, 
which is represented by an infinite sum over the partial 
waves. The numerical approach was developed in the pre- 
vious works (25U27| , where it was used for calculating the 
QED and electron-electron interaction corrections to the 
radiative recombination of electrons with highly charged 



For a given value of the relativistic angular momen- 
tum quantum number k, the radial part of the Dirac 
Coulomb Green function is represented in terms of the 
two-component solutions of the radial Dirac equation reg- 
ular at the origin (S®,) and the infinity (</>£°), 



G K (E, n, r 2 ) = - ^{E, n) 0f (E, r 2 ) 6{n - ra) 



(20) 



where E denotes the energy argument of the Green func- 
tion, r\ and r 2 are the radial arguments, and 9 is the step 
function. For the point Coulomb potential, the regular 
and irregular solutions of the radial Dirac equation are 
expressed analytically in terms of the Whittaker func- 
tions of the first and second kind (for explicit formulas 
see, e.g., Ref. (H). 

When the energy argument E is real and greater than 
the electron rest mass E > m, the Dirac Green function 
is a complex multi-valued function. So, care should be 
taken in this case to choose the correct branch (i.e., the 
sign of the imaginary part) of the Green function. The 
branch of the Green function is fixed by the sign of the 
infinitesimal imaginary addition iO in the energy denom- 
inators of Eq. (|16p . Specifically, the energy argument of 
the Green function in the first term in the brackets of 
Eq. (|16l) . £q a + oj, is greater than the electron mass. For 
positive-energy intermediate states, the energy denomi- 
nator takes the form eQ a +uj — e n +i0, which implies that 
the energy argument of the Green function has a small 
positive imaginary addition and, therefore, the branch 
cut of the Green function [m, 00) should be approached 
from above. 

A serious problem arises in the numerical evaluation of 
the radial integrals for the left graph in Fig. [1] (with the 
electron correlation modifying the final-state wave func- 
tion) . In this case, the continuum-state Dirac wave func- 
tion has to be integrated together with the Dirac Green 
function with the energy argument E > m and the spher- 
ical Bessel function. All three functions are strongly os- 
cillating and slowly decreasing for large radial arguments. 
It is practically impossible to accurately evaluate such in- 
tegral by a straightforward numerical integration. In this 
work, we use the method of the complex-plane rotation 
of the integration contour, which was previously applied 
by us to the evaluation of the free- free transition integrals 
in the bremsstrahlung (29[. 
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We now consider the evaluation of the problematic ra- 
dial integrals in more details. Let us introduce the radius 
of the atom 1Z, which is defined as the smallest distance 
at which all the bound-state electron wave functions van- 
ish. In the inner region r < 1Z, evaluation of each partial- 
wave expansion term of the amplitude [Eq. (|16[) ] involves 
a three-dimensional radial integration. In the outer re- 
gion r > 1Z, however, all the integrals involving the 
bound-state wave functions reach their asymptotical val- 
ues and only an one-dimensional integral of the free-free 
type needs to be evaluated. The general form of such 
integral is 

/>oo 

J K = / drr 2 j L (k 3 r)fl(E 2 , r) <f>£ (E 1: r) , (21) 

where jl is the spherical Bessel function originating from 
the photon propagator in the electron-electron interac- 
tion, f % is the radial component of the Dirac continuum 
wave function (i = 1,2), and (ft 003 is the radial component 
of the irregular solution of the Dirac equation originating 
from the Green function (J — 1,2). The energy conser- 
vation requires that E\ = E% + fca, which leads to the 

inequality p\ > p 2 + k%, where pi^ = yj E 2 2 — m 2 are 

the electron momenta associated with the corresponding 
energies. We now analytically continue the integrand of 
Eq. (|2ip into the complex r plane and take into account 
the asymptotical behaviour of individual functions for 
large values of Im(r), 

f K2 (E 2 ,r) ~exp[p 2 |Im(r)|] , (22) 
4>™(Ei,r) ~ exp [ Pl Im(r)] , (23) 
jj(fc 3 r) ~exp[fc 3 |Im(r)|] , (24) 

where only the leading exponential behaviour is kept. 
We observe that, rotating the integration contour in 
the integral (|2"Tj) into the lower complex half-plane r — > 
—ir, we transform the strongly oscillating integrand 
into an exponentially decreasing one, which falls off as 
exp [— (pi — P2 — k$) r] for large r. After the rotation of 
the contour, the integral can be easily evaluated numer- 
ically up to a desired precision. 

After the radial integrals are successfully evaluated, 
the remaining problem is the summation of the partial- 
wave expansions. Altogether there are five partial-wave 
expansions to be delt with: the expansions of the both fi- 
nal electron wave functions, the Green function, the wave 
function of the absorbed photon, and the photon propa- 
gator of the electron-electron interaction. After the an- 
gular momentum selection rules are taken into account, 
only two partial-wave expansions of the two final-state 
electrons remain unbound. The corresponding expansion 
parameters are K\ and k 2 in Eq. (1141) . The convergence 
of the resulting double partial-wave expansion is good 
for the photon energies near the threshold cj ~ cj cr and 
a non-symmetric energy sharing x <C V 2 but gradually 
deteriorates when oj increases and x approaches 1 /2. In 
the most difficult case considered here, w/w cr sa 30 and 



x = 1/2, up to (|«i|, |k 2 |) = (15, 15) partial waves were 
included into the calculation. The truncated partial wave 
expansion was extrapolated into infinity by using the e 
resummation algorithm (see, e.g., Ref. [30|). 

In order to check the numerical procedure, we eval- 
uated the nonrelativistic limit of our calculations. The 
easiest way to do this is to decrease the value of the 
fine structure constant (or, equivalently, to increase the 
value of the speed of light) by a large factor. However, 
a straightforward implementation of this scheme leads to 
numerical instabilities, since the relativistic operators in- 
volve products of the upper and the lower components of 
the Dirac wave function and the lower component van- 
ishes in the nonrelativistic limit. In order to avoid this 
problem, we express the lower component in terms of the 
upper component as 

/«W = ifx + — W)> ( 25 ) 
Zm \ar r J 

which is valid to the leading order in Za. After this sub- 
stitution, our relativistic code yielded a stable numerical 
limit when the fine structure constant was decreased by 
orders of magnitude, thus giving us the nonrelativistic 
limit. 



V. RESULTS AND DISCUSSION 

We start this section with a discussion of the total cross 
section of the double photoionization. As explained pre- 
viously, our approach accounts for the electron correla- 
tion to the first order of perturbation theory. Within 
this approximation we perform a rigorous relativistic 
treatment, without any further simplifications. This ap- 
proach is expected to yield accurate results for heavy 
helium-like ions, whereas for light ions, the precision of 
the method gradually deteriorates with decrease of the 
nuclear charge. For helium, the relativistic effects are 
weak but the electron correlation is strong, so that in 
this case our approach predictably yields worse results 
than the modern nonrelativistic methods. However, we 
will use the helium case for estimating the higher-order 
electron correlation effects omitted, having in mind that 
in helium-like ions these effects are suppressed by the 
inverse power of the nuclear charge. 

In Fig. [5J we compare the nonrelativistic limit of our 
calculations of the ratio of the double-to-single photoion- 
ization cross section R(u>) = a ++ /a + in helium-like ions 
with the experimental data available for helium and with 
the theoretical high-energy limits. In order to facili- 
tate the analysis of the results, we plot the scaled ratio 
Z 2 R(uj) as a function of the ratio oj/uj ci . 

Firstly, we confirm the known statement [3l|, |32| that 
the nonrelativistic limit of the scaled ratio R(u)) = 
Z 2 R{uj I uj cr ) is well described by a universal function that 
does not depend on the nuclear charge number. We will 
show later that this universal scaling is strongly violated 
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FIG. 2: (Color online) Comparison of the nonrelativistic limit 
of the present calculations with the experimental results for 
helium and the nonrelativistic theoretical high-energy lim- 
its. The solid (black) line shows our nonrelativistic results 
for Z — 2; the dashed (blue) line, for Z = 40; the dashed- 
dotted (red) line, for Z — 54. (The three lines are practically 
indistinguishable on the picture.) The triangles, squares, and 
diamonds represent the experimental results for helium ob- 
tained in Refs. 34-37}], respectively. The dotted line shows 
the nonrelativistic high-energy limit calculated to the leading 
order of the perturbation theory [3]. The dash-dot-dotted line 
shows the nonrelativistic asymptotical limit calculated with 
the fully correlated wave functions [!, [f| . 



by the relativistic effects. Fig. [2] shows that the numer- 
ical results calculated for different nuclear charges are 
practically indistinguishable from each other in the non- 
relativistic limit. 

Secondly, we observe good agreement between our non- 
relativistic results and the leading term of the 1/Z expan- 
sion of the asymptotics 0,111, R a/z = 0.0932/Z 2 . The 
deviation of our numerical results from the experimental 
data [34T - [37| is consistent with the deviation of the per- 
turbative asymptotical value Rouz from the fully corre- 
lated result Ro, corl - = 0.0658/Z 2 [5, 6]. We, therefore, es- 
timate the higher-order electron-correlation effects omit- 
ted in the present treatment to be about 2/Z x 30% for 
the ratio R(u>). 

We now turn to our relativistic calculations. Fig. [3] 
presents our numerical results obtained with a step-by- 
step inclusion of individual relativistic effects for several 
helium-like ions. The four different treatments compared 
are (i) the nonrelativistic calculation (dotted line), (ii) 
the calculation with the relativistic wave functions but 
with neglecting the retardation in the electron-electron 
interaction and all multipoles of the absorbed photon 
higher than the relativistic El transition (dash-dotted 
line), (iii) the calculation with the relativistic wave func- 
tions and the full retardation but without the higher mul- 
tipoles of the absorbed photon (dashed line) , and (iv) the 



full relativistic treatment (solid line). 

We observe that for moderate photon energies, 
cj/w cr < 3, the dominant relativistic effect is the retarda- 
tion in the electron-electron interaction, whereas for high 
photon energies, lo/oj ci > 10, the effect of the higher 
multipoles of the absorbed photon becomes dominant. 
Altogether, the relativistic effects are large and change 
qualitatively the shape of the energy dependence of the 
ratio R{<jj) already for medium- Z ions. In particular, the 
maximum of the curve located for helium at w/uj cr rj 2.5 
disappears for Z > 20 and R(oj) becomes a monotoni- 
cally growing function (at least, up to the maximal pho- 
ton energies accessible in our calculations), which is in 
contrast to the nonrelativistic case where R(oj) gradually 
decreases to approach a constant high-energy limit from 
above. It is interesting to note that similarly strong rel- 
ativistic effects were recently reported for the nonelastic 
electron scattering from the hydrogen- like ions [HJ . Sig- 
nificant effects caused by the magnetic dipole and electric 
quadrupole transitions in the double photoionization of 
helium-like ions were previously reported in Ref. [39| . 

In order to check our relativistic calculations, we em- 
ployed different gauges for the absorbed photon and, 
separately, for the photon propagator of the electron- 
electron interaction. We found that the gauge-dependent 
term in the photon absorbtion operator [see Eq. ([3])] 
vanishes in the actual calculations. Independently, we 
demonstrated that our calculations with the electron- 
electron interaction operator in the Feynman and the 
Coulomb gauges [see Eqs. (fl"3j) and ([!"§)) ] yield the same 
results. This was an important cross-check of the nu- 
merical procedure since the contributions of the two di- 
agrams in Fig. [T] are separately not gauge invariant. As 
an additional test, we calculated the contribution of the 
right graph in Fig. [1] independently by a different numeri- 
cal technique based on the -B-spline representation of the 
Dirac Coulomb Green function [4fj| . 

The results of our relativistic calculations for vari- 
ous helium-like ions are summarized in Fig. 3J It can 
be readily seen that the relativistic effects strongly vio- 
late the nonrelativistic scaling rule stating that the func- 
tion R(oj) = Z 2 R(uj/uj cr ) does not depend on the nu- 
clear charge. The relativistic enhancement of the ratio 
of the double-to-single photoionization cross sections can 
be conveniently parameterized as 



R{uj, Z) = i?NR,(w, Z) 



l + (Za) 2 frel(w,Z) 



(26) 



where i?NR is the nonrelativistic limit of the ratio 
(7 ++ /ct + and / re i is a smooth function of the nuclear 
charge Z and photon energy u. This function is plot- 
ted for several nuclear charges in Fig. [5] We observe 
that the numerical values of / re i lay within the interval 
of (1.5, 3.5) for the wide range of nuclear charge numbers 
and photon energies. We also found that for Z < 20, 
/ re l « 3 and is nearly Z-independent. In particular, for 
helium and 1 keV photon energy (w/u; cr « 13), the above 
formula predicts the relativistic enhancement of about 
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0.8%, which is negligible at the current level of experi- 
mental precision. 

The normalized energy distribution of the ejected elec- 
trons is shown in Fig. [5] for several nuclear charges and 
photon energies. Again, we observe that the scaled en- 
ergy distribution does not depend on the nuclear charge 
in the nonrelativistic limit and that this scaling is vio- 
lated by the relativistic effects. Note that in the case 
of the energy distribution, the differential cross section 
has to be scaled by the factor of Z 2 (Z a) 2 / a + and not by 
Z 2 /a + as for the total cross section; the additional factor 
is due to the fact that the integration interval in Eq. (TT5)) 
is proportional to (Za) 2 . We also confirm the known 
statement [I[ that the highly asymmetric energy sharing 
between the two emitted electrons dominates at the high 
photon energies, with one of the electrons carrying away 
most of the photon energy. 

VI. CONCLUSION 

In this paper, we presented a relativistic calculation of 
the double photoionization of helium-like atoms. Our ap- 
proach is based on the partial-wave representation of the 
Dirac continuum states and accounts for the retardation 
in the electron-electron interaction, the higher-order mul- 
tipoles of the absorbed photon as well as the interaction 
of the electrons with the nucleus without any expansion 
in the binding field. The electron-electron interaction 
is taken into account to the first order of perturbation 
theory. The omitted higher-order electron correlation ef- 
fects are estimated by comparing our numerical results 
for atomic helium with the experimental results and the 
available nonrelativistic theory. The calculational results 
are shown to be gauge invariant both with respect of 
the gauge of the absorbed photon and the gauge of the 
electron-electron interaction. 

Our calculation shows that the relativistic effects be- 
come prominent in the double photoionization cross sec- 
tion already for medium- Z ions. These effects change 
the shape of the energy dependence of the ratio of the 
double-to-single photoionization R(uj) drastically. In par- 
ticular, the well-known constant high-energy asymptotic 
behaviour of R(u>) in helium gives place to the monotoni- 



cally growing behaviour in the case of helium-like targets 
with Z > 20. 

We note that the reported magnitude of the relativis- 
tic effects is insufficient to explain the large discrepancy 
between theory and experiment observed for the double 
K-sheW photoionization in neutral atoms. In particular, 
for the photon energy of 90 keV in silver (ui/u; CI = 1.73), 
we obtain the relativistic enhancement of about 50% [see 
Eq. (|26p]. which is much smaller than the 4- fold disagree- 
ment between theory and experiment [13J. 

It should be mentioned that in the present investiga- 
tion we do not consider the Compton scattering mecha- 
nism of the total double photoionization, which is known 
to contribute significantly for high-energy photons [f|. 
This reaction channel is more difficult for a relativistic 
calculation than the photoabsorption considered in the 
present work. To a certain extent, the photoabsorption 
and the Compton mechanisms of the double photoioniza- 
tion can be distinguished in the experiment by measuring 
the momentum vector of the recoil ions fol ■ 
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Appendix A: Single photoionization: calculation 
formulas 

In order to perform integrations over the angular vari- 
ables in the matrix element of the photon absorption 
operator, we fix the z-axis of our coordinate system to 
be directed along the photon momentum k and assume 
that the polarization vector u has the only nonvanishing 
spherical component u\ (A = ±1). The result for the ma- 
trix element of the gauge-independent part of the photon 
absorption operator is 



( Ka ^ a \ a ■ u x e lk - r \ Kn Hn) = (-1) ]T i 1+L U L (-l)i--/^. C^ atjn _^ CH 1X P JL (an) , (Al) 

JL 

where II^ = \/2L + 1 and the initial and final states are the Dirac states with given relativistic angular momentum 
quantum number k and the momentum projection /i. The radial integrals Pjl are 

/>oo 

P JL (an)= \ drr 2 j L (uir) [g n (r) f a (r) S,j L (Kn,-K a ) - f,i(r) g a (r) S JL (~K n ,K a )] , (A2) 
Jo 

where is the spherical Bessel function, <?.j and fi are the upper and the lower radial components of the Dirac wave 
function, respectively, and Sjl are the angular coefficients given by Eqs. (C7)-(C9) of Ref. (42[. The matrix element 
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FIG. 3: (Color online) Comparison of the results of the fully relativistic calculations (solid line, black) with results of various 
approximate treatments. The dotted (red) line represents the nonrelativistic limit. The dashed-dotted (blue) line is obtained 
with the relativistic wave functions but without retardation and with including the relativistic El transition only. The dashed 
(green) line shows the results obtained with the relativistic wave functions and the full retardation but with including the 
relativistic El transition only. The threshold photon energy of the double photoionization is cj cr = 2.56 keV for Ne 8+ , 
uj cr = 10.60 keV for Ca 18+ , and oj cr = 43.75 keV for Zr 38+ . 
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FIG. 4: (Color online) Relativistic results for the scaled ratio 
of the double photoionization cross section a ++ to the sin- 
gle photoionization cross section a + , for helium-like ions with 
different nuclear charges Z, as a function of the photon en- 
ergy uj divided by the threshold photon energy of the double 
photoionization lj ct . 



of the gauge-dependent part of the photon absorption operator can be evaluated as 



(K a ^ a \(a -k-1) e ih - r \ Knf x n ) = (-1)^ L IT L (-l)^-Mn cH 



JL 



i C Lo,w p JL{an) + 5 JtL C L (n a , n n ) R L (an) , 

(A3) 
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FIG. 5: (Color online) Relativistic enhancement function / re i 
defined by Eq. (|26[) . for different helium- like ions. 
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FIG. 6: (Color online) Normalized energy distribution of the ejected electrons as a function of the energy sharing parameter x 
defined by Eqs. (|12[) and (|13|l . for different nuclear charge numbers and photon energies. 



where the angular coefficients Cl are given by Eq. (C10) of Ref. [42| and Rl is the radial integral 

/>oo 

R L {an)= / dr r 2 j L (wr) [g n (r ) g a (r ) + /„ (r ) f a (r)] . 
Jo 



(A4) 



r 



Appendix B: Double photoionization: calculation 
formulas 



Let us write the amplitude ([16]) as 

T A h+ = Ti2 Q & — T 2 i a ;, — Ti2ba 



1~21b 



(Bl) 



where Tyiab corresponds to the part of Eq. (|16p with 
Pe\Pei = an d QaQb — ab and the remaining 

three terms are obtained by permutations. The contri- 



butions due to the two terms in the braces of Eq. I|16p 
will be denoted by subscripts A and B, respectively, 

T 12ab — T 12ab,A + T 12ab,B ■ 

The calculation formulas take the simplest form in the 
case when the initial electron state is the ground state 
of the atom. In this case, the permutation over the a 
and b electrons yields just a combinatorial factor and 
all summations over the momentum projections can be 
evaluated in the closed form. The result for the single 
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differential cross section is 
da ++ An 2 a 



The amplitude is r 12ab = T 12a b,A + T 12a b,B, 



etei 



(B2) 



n 2aM - ^ (-i) , n L c i0 , 1A | i . n L , 1 ^ — 



c a + cj - e n 



(B3) 



LL'. 

and 

a i ix,i+Ltt /^ja s j n ,j2 PjL(ein)R L '(A,ne 2 ab) 



Tl 2a b,B 

LL 



4= y (-i) ^ u y rjLw>"VK~,>^™> (B4) 



r 



where i?i(A, abed) is the relativistic generalization of the by Eqs. (C1)-(C10) of Ref. [42[ and A = e 2 — £&. 
Slater integral for the electron-electron interaction given 
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